The PDE4 Inhibitors Roflumilast and Rolipram Rescue ADO2 Osteoclast Resorption Dysfunction.

Autosomal Dominant Osteopetrosis type II (ADO2) is a rare bone disease of impaired osteoclastic bone resorption caused by heterozygous missense mutations in the chloride channel 7 (CLCN7). Adenylate cyclase, which catalyzes the formation of cAMP, is critical for lysosomal acidification in osteoclasts. We found reduced cAMP levels in ADO2 osteoclasts compared to wild-type (WT) osteoclasts, leading us to examine whether regulating cAMP would improve ADO2 osteoclast activity. Although forskolin, a known activator of adenylate cyclase and cAMP levels, negatively affected osteoclast number, it led to an overall increase in ADO2 and WT osteoclast resorption activity in vitro. Next, we examined cAMP hydrolysis by the phosphodiesterase 4 (PDE4) proteins in ADO2 versus WT osteoclasts. QPCR analysis revealed higher expression of the three major PDE4 subtypes (4a, 4b, 4d) in ADO2 osteoclasts compared in WT, consistent with reduced cAMP levels in ADO2 osteoclasts. In addition, we found that the PDE4 antagonists, rolipram and roflumilast, stimulated ADO2 and WT osteoclast formation in a dose-dependent manner. Importantly, roflumilast and rolipram displayed a concentration-dependent increase in osteoclast resorption activity which was greater in ADO2 than WT osteoclasts. Moreover, treatment with roflumilast rescued cAMP levels in ADO2 OCLs. The key findings from our studies demonstrate that osteoclasts from ADO2 mice exhibit reduced cAMP levels and PDE4 inhibition rescues cAMP levels and ADO2 osteoclast activity dysfunction in vitro. The mechanism of action of PDE4 inhibitors and their ability to reduce the high bone mass of ADO2 mice in vivo are currently under investigation. Importantly, these studies advance the understanding of the mechanisms underlying the ADO2 osteoclast dysfunction which is critical for the development of therapeutic approaches to treat clinically affected ADO2 patients.

Effect of Roflumilast, a Selective PDE4 Inhibitor, on Bone Phenotypes in ADO2 Mice.

Autosomal Dominant Osteopetrosis type II (ADO2) is a rare bone disease of impaired osteoclastic bone resorption that usually results from heterozygous missense mutations in the chloride channel 7 (CLCN7) gene. We previously created mouse models of ADO2 (p.G213R) with one of the most common mutations (G215R) as found in humans and demonstrated that this mutation in mice phenocopies the human disease of ADO2. Previous studies have shown that roflumilast (RF), a selective phosphodiesterase 4 (PDE4) inhibitor that regulates the cAMP pathway, can increase osteoclast activity. We also observed that RF increased bone resorption in both wild-type and ADO2 heterozygous osteoclasts in vitro, suggesting it might rescue bone phenotypes in ADO2 mice. To test this hypothesis, we administered RF-treated diets (0, 20 and 100 mg/kg) to 8-week-old ADO2 mice for 6 months. We evaluated bone mineral density and bone micro-architecture using longitudinal in-vivo DXA and micro-CT at baseline, and 6-, 12-, 18-, and 24-week post-baseline time points. Additionally, we analyzed serum bone biomarkers (CTX, TRAP, and P1NP) at baseline, 12-, and 24-week post-baseline. Our findings revealed that RF treatment did not improve aBMD (whole body, femur, and spine) and trabecular BV/TV (distal femur) in ADO2 mice compared to the control group treated with a normal diet. Furthermore, we did not observe any significant changes in serum levels of bone biomarkers due to RF treatment in these mice. Overall, our results indicate that RF does not rescue the osteopetrotic bone phenotypes in ADO2 heterozygous mice.

Autosomal Dominant Osteopetrosis (ADO) caused by a Missense Variant in the TCIRG1 Gene.

CONTEXT: Autosomal dominant osteopetrosis (ADO) is a rare genetic disorder due to impaired osteoclastic bone resorption. Clinical manifestations frequently include fractures, osteonecrosis (particularly of the jaw or maxilla), osteomyelitis, blindness, and/or bone marrow failure. ADO usually results from heterozygous missense variants in the Chloride Channel 7 gene (CLCN7) that cause disease by a dominant negative mechanism. Variants in the T cell immune regulator 1 gene (TCIRG1) are commonly identified in autosomal recessive osteopetrosis but have only been reported in one patient with ADO. CASE DESCRIPTION: Here we report 3 family members with a single heterozygous missense variant (p.Gly579Arg) in TCIRG1 who have a phenotype consistent with ADO. Three of five protein prediction programs suggest this variant likely inhibits the function of TCIRG1. CONCLUSIONS: This is the first description of adult presentation of ADO caused by a TCIRG1 variant. Similar to families with ADO from CLCN7 mutations, this variant in TCIRG1 results in marked phenotype variability, with two subjects having severe disease and the third having very mild disease. This family report implicates TCIRG1 missense mutations as a cause of ADO and demonstrates that the marked phenotypic variability in ADO may extend to disease caused by TCIRG1 missense mutations.

Bone marrow transplantation as a therapy for autosomal dominant osteopetrosis type 2 in mice.

Autosomal dominant osteopetrosis type II (ADO2) is a heritable bone disease of impaired osteoclastic bone resorption caused by missense mutations in the chloride channel 7 (CLCN7) gene. Clinical features of ADO2 include fractures, osteomyelitis of jaw, vision loss, and in severe cases, bone marrow failure. Currently, there is no effective therapy for ADO2, and patients usually receive symptomatic treatments. Theoretically, bone marrow transplantation (BMT), which is commonly used in recessive osteopetrosis, could be used to treat ADO2, although the frequency of complications related to BMT is quite high. We created an ADO2 knock-in (p.G213R mutation) mouse model on the 129 genetic background, and their phenotypes mimic the human disease of ADO2. To test whether BMT could restore osteoclast function and rescue the bone phenotypes in ADO2 mice, we transplanted bone marrow cells from 6-8 weeks old male WT donor mice into recipient female ADO2 mice. Also, to determine whether age at the time of transplant may play a role in transplant success, we performed BMT in young (12-week-old) and old (9-month-old) ADO2 mice. Our data indicate that ADO2 mice transplanted with WT marrow achieved more than 90% engraftment up to 6 months post-transplantation at both young and old ages. The in-vivo DXA data revealed that young ADO2 mice transplanted with WT marrow had significantly lower whole body and spine areal bone mineral density (aBMD) at month 6 post-transplantation compared to the ADO2 control mice. The old ADO2 mice also displayed significantly lower whole body, femur, and spine aBMD at months 4 and 5 post-transplantation compared to the age-matched control mice. The in-vivo micro-CT data showed that ADO2 experimental mice transplanted with WT marrow had significantly lower BV/TV at months 2 and 4 post-transplantation compared to the ADO2 control mice at a young age. In contrast, ADO2 control and experimental mice displayed similar BV/TV values for all post-transplantation time points at old age. In addition, serum CTX was significantly higher at month 2 post-transplantation in both young and old ADO2 experimental mice compared to the ADO2 control mice. Serum P1NP levels in young ADO2 experimental mice were significantly higher at baseline and month 2 post-transplantation compared to the ADO2 control mice. These data suggest that BMT may provide, at least, some beneficial effect at both young and adult ages.

Unusual Cortical Phenotype After Hematopoietic Stem Cell Transplantation in a Patient With Osteopetrosis.

The osteopetroses are a group of rare genetic diseases caused by osteoclast dysfunction or absence. The hallmark of osteopetrosis is generalized increased bone mineral density (BMD). However, the bone is fragile and fractures are common. Autosomal recessive osteopetrosis is usually a severe disorder and often life-threatening in childhood. We present male siblings with autosomal recessive osteopetrosis due to biallelic variants in TCIRG1 who survived childhood and underwent hematopoietic stem cell transplant (HSCT) in adulthood. One sibling died of posttransplant complications. After transplant, the other sibling had improvement of multiple clinical parameters, including some decline in BMD Z-scores by dual-energy X-ray absorptiometry (DXA) and cessation of fractures. However, spine quantitative computed tomography 11 years after transplant demonstrated an anvil pattern of sclerosis with BMD Z-score of +18.3. High-resolution peripheral quantitative computed tomography (HR-pQCT) of the tibia demonstrated near complete obliteration of the marrow space combined with an unusual cortical phenotype, suggesting extensive cortical porosity at the distal tibia. This case highlights that despite successful transplantation and subsequent improvement in clinical parameters, this patient continued to have significantly elevated bone density and decreased marrow space. Transplant-associated increased cortical porosity is multifactorial and occurs in two-thirds of non-osteopetrotic patients undergoing HSCT. This finding after transplant in osteopetrosis may suggest particular sensitivity of the cortical bone to resorptive activity of transplanted osteoclasts. The case also suggests HR-pQCT may be a useful modality for imaging and assessing the therapeutic effects on bone in individuals with osteopetrosis. © 2022 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Chloroquine increases osteoclast activity in vitro but does not improve the osteopetrotic bone phenotype of ADO2 mice.

Autosomal Dominant Osteopetrosis type II (ADO2) is a bone disease of impaired osteoclastic bone resorption that usually results from heterozygous missense mutations in the chloride channel 7 (CLCN7) gene. We created mouse models of ADO2 by introducing a knock-in (p.G213R) mutation in the Clcn7 gene, which is analogous to one of the common mutations (G215R) found in humans. The mutation leads to severe osteopetrosis and lethality in homozygous mice but produces substantial phenotypic variability in heterozygous mice on different genetic backgrounds that phenocopy the human disease of ADO2. ADO2 is an osteoclast-intrinsic disease, and lysosomal enzymes and proteins are critical for osteoclast activity. Chloroquine (CQ) is known to affect lysosomal trafficking, intracellular signaling and the lysosomal and vesicular pH, suggesting it might improve ADO2 osteoclast function. We tested this hypothesis in cell culture studies using osteoclasts derived from wild-type (WT or ADO2+/+) and ADO2 heterozygous (ADO2+/-) mice and found that CQ and its metabolite desethylchloroquine (DCQ), significantly increased ADO2+/- osteoclasts bone resorption activity in vitro, whereas bone resorption of ADO2+/+ osteoclasts was increased only by DCQ. In addition, we exploited our unique animal model of ADO2 on 129 background to identify the effect of CQ for the treatment of ADO2. Female ADO2 mice at 8 weeks of age were treated with 5 doses of CQ (1, 2.5, 5, 7.5 and 10 mg/kg BW/day) via drinking water for 6 months. Bone mineral density and bone micro-architecture were analyzed by longitudinal in vivo DXA and micro-CT at baseline, 3 and 6 months. Serum bone biomarkers (CTX, TRAP and P1NP) were also analyzed at these time points. CQ treatment at the doses tested failed to produce any significant changes of aBMD, BMC (whole body, femur and spine) and trabecular BV/TV (distal femur) in ADO2 mice compared to the control group (water only). Further, levels of bone biomarkers were not significantly changed due to CQ treatment in these mice. Our findings indicate that while CQ increased osteoclast activity in vitro, it did not improve the osteopetrotic bone phenotypes in ADO2 heterozygous mice.

Overexpression of WNT16 Does Not Prevent Cortical Bone Loss Due to Glucocorticoid Treatment in Mice.

Glucocorticoids (GC) are commonly used for the treatment of a wide variety of autoimmune, pulmonary, gastrointestinal, and malignancy conditions. One of the devastating side effects of GC use is osteoporotic fractures, particularly in the spine and hip. Bisphosphonates (BP) are the most commonly prescribed pharmacological agents for the prevention and treatment of GC-induced osteoporosis (GIO). However, GIO is marked by reduced bone formation and BP serves mainly to decrease bone resorption. The WNT signaling pathway plays a major role in bone and mineral homeostasis. Previously, we demonstrated that overexpression of WNT16 in mice led to higher bone mineral density and improved bone microarchitecture and strength. We hypothesized that WNT16 overexpression would prevent bone loss due to glucocorticoid treatment in mice. To test our hypothesis, we treated adult wild-type and WNT16-transgenic mice with vehicle and GC (prednisolone; 2.1 mg/kg body weight) via slow-release pellets for 28 days. We measured bone mass and microarchitecture by dual-energy X-ray absorptiometry (DXA) and micro-CT, and performed gene expression and serum biochemical analysis. We found that GC treatment compared with the vehicle significantly decreased femoral areal bone mineral density (aBMD), bone mineral content (BMC), and cortical bone area and thickness in both wild-type and transgenic female mice. In contrast, the trabecular bone parameters at distal femur were not significantly changed by GC treatment in male and female mice for both genotypes. Further, we observed significantly lower level of serum P1NP and a tendency of higher level of serum TRAP in wild-type and transgenic mice due to GC treatment in both sexes. Gene expression analysis showed lower mRNA levels of Wnt16, Opg, and Opg/Rankl ratio in GC-treated female mice for both genotypes compared with the sex-matched vehicle-treated mice. These data suggest that although WNT16 overexpression resulted in higher baseline bone mineral density and bone volume per trabecular volume (BV/TV) in the transgenic mice, this was insufficient to prevent bone loss in mice due to glucocorticoid treatment.

Phenotypic severity of autosomal dominant osteopetrosis type II (ADO2) mice on different genetic backgrounds recapitulates the features of human disease.

Autosomal dominant osteopetrosis type II (ADO2) is a heritable osteosclerotic bone disorder due to dysfunctional osteoclast activity. ADO2 is caused by missense mutations in the chloride channel 7 (CLCN7) gene characterized by osteosclerosis with multiple fractures. ADO2 can result in osteomyelitis, visual loss and bone marrow failure. Currently, there is no cure for ADO2, and until recently no appropriate animal model of ADO2 existed to understand better the pathogenesis of this disease and to test new therapies. Therefore, we created ADO2 knock-in mouse model with a G213R (human homolog of G215R) missense mutation in the Clcn7 gene on 129S1 background, and demonstrated that this mouse model phenocopies human ADO2. As ADO2 gives rise to incomplete penetrance (66%) in human and marked phenotypic variability is observed among patients with the same mutation, we hypothesized that the severity and penetrance of ADO2 will also vary in mouse models on different genetic backgrounds. To test this, we created ADO2 mouse models in DBA/D2, C57BL/6J/B6 and Balb/c strains, and compared bone phenotypes and performed serum biochemical analysis between strain- and age-matched wild-type (WT) and ADO2 mice. At 3months of age, whole body aBMD was higher (4-7% in male; 1-5% in female) in the ADO2 mice compared to their wild-type littermates. In addition, ADO2 male mice on 129 background displayed highest percent increase of BV/TV (106%), followed by D2 (92%), B6 (46%), and Balb/c (33%) compared to strain-matched wild-type mice. We observed similar differences for BV/TV between ADO2 and wild-type mice on different genetic backgrounds in female: 129 (96%)>D2 (73%)>Balb/c (39%) and B6 (36%). Serum calcium, phosphorus, alkaline phosphatase and P1NP levels were similar in the WT and ADO2 mice on all genetic backgrounds but TRAP was higher (76% to 220% in male; 33-95% in female) and CTX/TRAP ratio was lower (39-65% in male and 3-41% in female) in the ADO2 mice compared to their strain-matched wild-type littermates. We also found that young (3months) ADO2 mice on 129S1 background exhibited 200% higher trabecular BV/TV whereas old (18months) ADO2 mice displayed 400-700% higher BV/TV compared to their age-matched wild-type controls. In summary, phenotypic severity in ADO2 mice varied markedly on different genetic backgrounds (129>D2>Balb/c>B6) and became more pronounced with age, which resembles the wide variations in phenotype observed in ADO2 patients. These mouse models will help us to identify genes/factors that influence severity and penetrance of ADO2, and test innovative therapies to treat this disease.

Bone Mass and Strength are Significantly Improved in Mice Overexpressing Human WNT16 in Osteocytes.

Recently, we demonstrated that osteoblast-specific overexpression of human WNT16 increased both cortical and trabecular bone mass and structure in mice. To further identify the cell-specific role of Wnt16 in bone homeostasis, we created transgenic (TG) mice overexpressing human WNT16 in osteocytes using Dmp1 promoter (Dmp1-hWNT16 TG) on C57BL/6 (B6) background. We analyzed bone phenotypes and serum bone biomarkers, performed gene expression analysis and measured dynamic bone histomorphometry in Dmp1-hWNT16 TG and wild-type (WT) mice. Compared to WT mice, Dmp1-hWNT16 TG mice exhibited significantly higher whole-body, spine and femoral aBMD, BMC and trabecular (BV/TV, Tb.N, and Tb.Th) and cortical (bone area and thickness) parameters in both male and female at 12 weeks of age. Femur stiffness and ultimate force were also significantly improved in the Dmp1-hWNT16 TG female mice, compared to sex-matched WT littermates. In addition, female Dmp1-hWNT16 TG mice displayed significantly higher MS/BS, MAR and BFR/BS compared to the WT mice. Gene expression analysis demonstrated significantly higher mRNA level of Alp in both male and female Dmp1-hWNT16 TG mice and significantly higher levels of Osteocalcin, Opg and Rankl in the male Dmp1-hWNT16 TG mice in bone tissue compared to sex-matched WT mice. These results indicate that WNT16 plays a critical role for acquisition of both cortical and trabecular bone mass and strength. Strategies designed to use WNT16 as a target for therapeutic interventions will be valuable to treat osteoporosis and other low bone mass conditions.

Osteoblast-Specific Overexpression of Human WNT16 Increases Both Cortical and Trabecular Bone Mass and Structure in Mice.

Previous genome-wide association studies have identified common variants in genes associated with bone mineral density (BMD) and risk of fracture. Recently, we identified single nucleotide polymorphisms (SNPs) in Wingless-type mouse mammary tumor virus integration site (WNT)16 that were associated with peak BMD in premenopausal women. To further identify the role of Wnt16 in bone mass regulation, we created transgenic (TG) mice overexpressing human WNT16 in osteoblasts. We compared bone phenotypes, serum biochemistry, gene expression, and dynamic bone histomorphometry between TG and wild-type (WT) mice. Compared with WT mice, WNT16-TG mice exhibited significantly higher whole-body areal BMD and bone mineral content (BMC) at 6 and 12 weeks of age in both male and female. Microcomputer tomography analysis of trabecular bone at distal femur revealed 3-fold (male) and 14-fold (female) higher bone volume/tissue volume (BV/TV), and significantly higher trabecular number and trabecular thickness but lower trabecular separation in TG mice compared with WT littermates in both sexes. The cortical bone at femur midshaft also displayed significantly greater bone area/total area and cortical thickness in the TG mice in both sexes. Serum biochemistry analysis showed that male TG mice had higher serum alkaline phosphatase, osteocalcin, osteoprotegerin (OPG), OPG to receptor activator of NF-kB ligand (tumor necrosis family ligand superfamily, number 11; RANKL) ratio as compared with WT mice. Also, lower carboxy-terminal collagen cross-link (CTX) to tartrate-resistant acid phosphatase 5, isoform b (TRAPc5b) ratio was observed in TG mice compared with WT littermates in both male and female. Histomorphometry data demonstrated that both male and female TG mice had significantly higher cortical and trabecular mineralizing surface/bone surface and bone formation rate compared with sex-matched WT mice. Gene expression analysis demonstrated higher expression of Alp, OC, Opg, and Opg to Rankl ratio in bone tissue in the TG mice compared with WT littermates. Our data indicate that WNT16 is critical for positive regulation of both cortical and trabecular bone mass and structure and that this molecule might be targeted for therapeutic interventions to treat osteoporosis.

Fine mapping of bone structure and strength QTLs in heterogeneous stock rat.

We previously demonstrated that skeletal structure and strength phenotypes vary considerably in heterogeneous stock (HS) rats. These phenotypes were found to be strongly heritable, suggesting that the HS rat model represents a unique genetic resource for dissecting the complex genetic etiology underlying bone fragility. The purpose of this study was to identify and localize genes associated with bone structure and strength phenotypes using 1524 adult male and female HS rats between 17 to 20 weeks of age. Structure measures included femur length, neck width, head width; femur and lumbar spine (L3-5) areas obtained by DXA; and cross-sectional areas (CSA) at the midshaft, distal femur and femoral neck, and the 5th lumbar vertebra measured by CT. In addition, measures of strength of the whole femur and femoral neck were obtained. Approximately 70,000 polymorphic SNPs distributed throughout the rat genome were selected for genotyping, with a mean linkage disequilibrium coefficient between neighboring SNPs of 0.95. Haplotypes were estimated across the entire genome for each rat using a multipoint haplotype reconstruction method, which calculates the probability of descent at each locus from each of the 8 HS founder strains. The haplotypes were then tested for association with each structure and strength phenotype via a mixed model with covariate adjustment. We identified quantitative trait loci (QTLs) for structure phenotypes on chromosomes 3, 8, 10, 12, 17 and 20, and QTLs for strength phenotypes on chromosomes 5, 10 and 11 that met a conservative genome-wide empiric significance threshold (FDR=5%; P<3×10(-6)). Importantly, most QTLs were localized to very narrow genomic regions (as small as 0.3 Mb and up to 3 Mb), each harboring a small set of candidate genes, both novel and previously shown to have roles in skeletal development and homeostasis.

Interferon Gamma, but not Calcitriol Improves the Osteopetrotic Phenotypes in ADO2 Mice.

ADO2 is a heritable osteosclerotic disorder that usually results from heterozygous missense dominant negative mutations in the chloride channel 7 gene (CLCN7). ADO2 is characterized by a wide range of features and severity, including multiple fractures, impaired vision due to secondary bony overgrowth and/or the lack of the optical canal enlargement with growth, and osteonecrosis/osteomyelitis. The disease is presently incurable, although anecdotal evidence suggests that calcitriol and interferon gamma-1b (IFN-G) may have some beneficial effects. To identify the role of these drugs for the treatment of ADO2, we utilized a knock-in (G213R mutation in Clcn7) ADO2 mouse model that resembles the human disease. Six-week-old ADO2 heterozygous mice were administered vehicle (PBS) or calcitriol or IFN-G 5 times per week for 8 weeks. We determined bone phenotypes using DXA and µCT, and analyzed serum biochemistry and bone resorption markers. ADO2 mice treated with all doses of IFN-G significantly (p<0.05) attenuated the increase of whole body aBMD and distal femur BV/TV gain in both male and female compared to the vehicle group. In contrast, mice treated with low and medium doses of calcitriol showed a trend of higher aBMD and BV/TV whereas high dose calcitriol significantly (p<0.05) increased bone mass compared to the vehicle group. The calcium and phosphorus levels did not differ between vehicle and IFN-G or calcitriol treated mice; however, we detected significantly (p<0.05) elevated levels of CTX/TRAP5b ratio in IFN-G treated mice. Our findings indicate that while IFN-G at all doses substantially improved the osteopetrotic phenotypes in ADO2 heterozygous mice, calcitriol treatment at any dose did not improve the phenotype and at high dose further increased bone mass. Thus, use of high dose calcitriol therapy in ADO2 patients merits serious reconsideration. Importantly, our data support the prospect of a clinical trial of IFN-G in ADO2 patients.

SIBLING family genes and bone mineral density: association and allele-specific expression in humans.

Osteoporosis is a common complex disorder with reduced bone mineral density (BMD) and increased susceptibility to fracture. Peak BMD is one of the primary determinants of osteoporotic fracture risk, and is under substantial genetic control. Extracellular matrix, a major component of the bone, influences BMD by regulating mineral deposition and maintaining cellular activity. It contains several SIBLING family proteins, null mutations of which cause mineralization defects in humans. In this study, we tested 59 single-nucleotide polymorphisms (SNPs) located in the 5 SIBLING family genes (DSPP, DMP1, IBSP, MEPE and SPP1) for association with normal variation in peak BMD in healthy men and women. We measured femoral neck (FN) and lumbar spine (LS) areal BMD by dual energy x-ray absorptiometry (DXA) in 1692 premenopausal European-American women, 512 premenopausal African-American women and 715 European-American men. SNPs were tested for association with FN and LS-BMD in the 3 subsamples. In the European-American women, we observed association (p≤0.005) with LS-BMD for SNPs in DSPP, IBSP and MEPE, and for FN-BMD with SNPs in DMP1 and IBSP. Allele-specific regulation of gene expression (ASE) is an important mechanism in which an allele giving rise to modest influence in transcript abundance might result in a predisposition to disease. To identify whether there was ASE of SIBLING family genes at these SNPs, we examined 52 human bone samples obtained from the femoral neck during surgical hip replacement (27 female, 25 male; 44 European-American and 8 African-American). We observed unidirectional ASE for the IBSP gene, with lower expression of the G allele compared to the A allele for SNP rs17013181. Our data suggest that SNPs within the SIBLING genes may contribute to normal variation of peak BMD. Further studies are necessary to identify the functional variants and to determine the mechanisms underlying the differences in ASE and how these differences relate to the pathophysiology of osteoporosis.

High-resolution genome screen for bone mineral density in heterogeneous stock rat.

We previously demonstrated that skeletal mass, structure, and biomechanical properties vary considerably in heterogeneous stock (HS) rat strains. In addition, we observed strong heritability for several of these skeletal phenotypes in the HS rat model, suggesting that it represents a unique genetic resource for dissecting the complex genetics underlying bone fragility. The purpose of this study was to identify and localize genes associated with bone mineral density in HS rats. We measured bone phenotypes from 1524 adult male and female HS rats between 17 and 20 weeks of age. Phenotypes included dual-energy X-ray absorptiometry (DXA) measurements for bone mineral content and areal bone mineral density (aBMD) for femur and lumbar spine (L3-L5), and volumetric BMD measurements by CT for the midshaft and distal femur, femur neck, and fifth lumbar vertebra (L5). A total of 70,000 polymorphic single-nucleotide polymorphisms (SNPs) distributed throughout the genome were selected from genotypes obtained from the Affymetrix rat custom SNPs array for the HS rat population. These SNPs spanned the HS rat genome with a mean linkage disequilibrium coefficient between neighboring SNPs of 0.95. Haplotypes were estimated across the entire genome for each rat using a multipoint haplotype reconstruction method, which calculates the probability of descent for each genotyped locus from each of the eight founder HS strains. The haplotypes were tested for association with each bone density phenotype via a mixed model with covariate adjustment. We identified quantitative trait loci (QTLs) for BMD phenotypes on chromosomes 2, 9, 10, and 13 meeting a conservative genomewide empiric significance threshold (false discovery rate [FDR] = 5%; p < 3 × 10(-6)). Importantly, most QTLs were localized to very small genomic regions (1-3 megabases [Mb]), allowing us to identify a narrow set of potential candidate genes including both novel genes and genes previously shown to have roles in skeletal development and homeostasis.

Generation of the first autosomal dominant osteopetrosis type II (ADO2) disease models.

Autosomal dominant osteopetrosis type II (ADO2) is a heritable osteosclerotic disorder dependent on osteoclast impairment. In most patients it results from heterozygous missense mutations in the chloride channel 7 (CLCN7) gene, encoding for a 2Cl(-)/1H(+) antiporter. By a knock-in strategy inserting a missense mutation in the Clcn7 gene, our two research groups independently generated mouse models of ADO2 on different genetic backgrounds carrying the homolog of the most frequent heterozygous mutation (p.G213R) in the Clcn7 gene found in humans. Our results demonstrate that the heterozygous model holds true presenting with higher bone mass, increased numbers of poorly resorbing osteoclasts and a lethal phenotype in the homozygous state. Considerable variability is observed in the heterozygous mice according with the mouse background, suggesting that modifier genes could influence the penetrance of the disease gene.

Combined sequence-based and genetic mapping analysis of complex traits in outbred rats.

Genetic mapping on fully sequenced individuals is transforming understanding of the relationship between molecular variation and variation in complex traits. Here we report a combined sequence and genetic mapping analysis in outbred rats that maps 355 quantitative trait loci for 122 phenotypes. We identify 35 causal genes involved in 31 phenotypes, implicating new genes in models of anxiety, heart disease and multiple sclerosis. The relationship between sequence and genetic variation is unexpectedly complex: at approximately 40% of quantitative trait loci, a single sequence variant cannot account for the phenotypic effect. Using comparable sequence and mapping data from mice, we show that the extent and spatial pattern of variation in inbred rats differ substantially from those of inbred mice and that the genetic variants in orthologous genes rarely contribute to the same phenotype in both species.

Heterogeneous stock rat: a unique animal model for mapping genes influencing bone fragility.

Previously, we demonstrated that skeletal mass, structure and biomechanical properties vary considerably among 11 different inbred rat strains. Subsequently, we performed quantitative trait loci (QTL) analysis in four inbred rat strains (F344, LEW, COP and DA) for different bone phenotypes and identified several candidate genes influencing various bone traits. The standard approach to narrowing QTL intervals down to a few candidate genes typically employs the generation of congenic lines, which is time consuming and often not successful. A potential alternative approach is to use a highly genetically informative animal model resource capable of delivering very high resolution gene mapping such as Heterogeneous stock (HS) rat. HS rat was derived from eight inbred progenitors: ACI/N, BN/SsN, BUF/N, F344/N, M520/N, MR/N, WKY/N and WN/N. The genetic recombination pattern generated across 50 generations in these rats has been shown to deliver ultra-high even gene-level resolution for complex genetic studies. The purpose of this study is to investigate the usefulness of the HS rat model for fine mapping and identification of genes underlying bone fragility phenotypes. We compared bone geometry, density and strength phenotypes at multiple skeletal sites in HS rats with those obtained from five of the eight progenitor inbred strains. In addition, we estimated the heritability for different bone phenotypes in these rats and employed principal component analysis to explore relationships among bone phenotypes in the HS rats. Our study demonstrates that significant variability exists for different skeletal phenotypes in HS rats compared with their inbred progenitors. In addition, we estimated high heritability for several bone phenotypes and biologically interpretable factors explaining significant overall variability, suggesting that the HS rat model could be a unique genetic resource for rapid and efficient discovery of the genetic determinants of bone fragility.

Identification of genes influencing skeletal phenotypes in congenic P/NP rats.

We previously showed that alcohol-preferring (P) rats have higher bone density than alcohol-nonpreferring (NP) rats. Genetic mapping in P and NP rats identified a major quantitative trait locus (QTL) between 4q22 and 4q34 for alcohol preference. At the same location, several QTLs linked to bone density and structure were detected in Fischer 344 (F344) and Lewis (LEW) rats, suggesting that bone mass and strength genes might cosegregate with genes that regulate alcohol preference. The aim of this study was to identify the genes segregating for skeletal phenotypes in congenic P and NP rats. Transfer of the NP chromosome 4 QTL into the P background (P.NP) significantly decreased areal bone mineral density (aBMD) and volumetric bone mineral density (vBMD) at several skeletal sites, whereas transfer of the P chromosome 4 QTL into the NP background (NP.P) significantly increased bone mineral content (BMC) and aBMD in the same skeletal sites. Microarray analysis from the femurs using Affymetrix Rat Genome arrays revealed 53 genes that were differentially expressed among the rat strains with a false discovery rate (FDR) of less than 10%. Nine candidate genes were found to be strongly correlated (r(2) > 0.50) with bone mass at multiple skeletal sites. The top three candidate genes, neuropeptide Y (Npy), alpha synuclein (Snca), and sepiapterin reductase (Spr), were confirmed using real-time quantitative PCR (qPCR). Ingenuity pathway analysis revealed relationships among the candidate genes related to bone metabolism involving beta-estradiol, interferon-gamma, and a voltage-gated calcium channel. We identified several candidate genes, including some novel genes on chromosome 4 segregating for skeletal phenotypes in reciprocal congenic P and NP rats.

Genes influencing spinal bone mineral density in inbred F344, LEW, COP, and DA rats.

Previously, we identified the regions of chromosomes 10q12-q31 and 15p16-q21 harbor quantitative trait loci (QTLs) for lumbar volumetric bone mineral density (vBMD) in female F2 rats derived from Fischer 344 (F344) x Lewis (LEW) and Copenhagen 2331 (COP) x Dark Agouti (DA) crosses. The purpose of this study is to identify the candidate genes within these QTL regions contributing to the variation in lumbar vBMD. RNA was extracted from bone tissue of F344, LEW, COP, and DA rats. Microarray analysis was performed using Affymetrix Rat Genome 230 2.0 Arrays. Genes differentially expressed among the rat strains were then ranked based on the strength of the correlation with lumbar vBMD in F2 animals derived from these rats. Quantitative PCR (qPCR) analysis was performed to confirm the prioritized candidate genes. A total of 285 genes were differentially expressed among all strains of rats with a false discovery rate less than 10%. Among these genes, 18 candidate genes were prioritized based on their strong correlation (r (2) > 0.90) with lumbar vBMD. Of these, 14 genes (Akap1, Asgr2, Esd, Fam101b, Irf1, Lcp1, Ltc4s, Mdp-1, Pdhb, Plxdc1, Rabep1, Rhot1, Slc2a4, Xpo4) were confirmed by qPCR. We identified several novel candidate genes influencing spinal vBMD in rats.

Differentially expressed genes strongly correlated with femur strength in rats.

The region of chromosome 1q33-q54 harbors quantitative trait loci (QTL) for femur strength in COPxDA and F344xLEW F2 rats. The purpose of this study is to identify the genes within this QTL region that contribute to the variation in femur strength. Microarray analysis was performed using RNA extracted from femurs of COP, DA, F344 and LEW rats. Genes differentially expressed in the 1q33-q54 region among these rat strains were then ranked based on the strength of correlation with femur strength in F2 animals derived from these rats. A total of 214 genes in this QTL region were differentially expressed among all rat strains, and 81 genes were found to be strongly correlated (r(2)>0.50) with femur strength. Of these, 12 candidate genes were prioritized for further validation, and 8 of these genes (Ifit3, Ppp2r5b, Irf7, Mpeg1, Bloc1s2, Pycard, Sec23ip, and Hps6) were confirmed by quantitative PCR (qPCR). Ingenuity Pathway Analysis suggested that these genes were involved in interferon alpha, nuclear factor-kappa B (NFkB), extracellular signal-related kinase (ERK), hepatocyte nuclear factor 4 alpha (HNF4A) and tumor necrosis factor (TNF) pathways.